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ABSTRACT 


Characteristics of broadband ELF, VLF, and LF emissions in the 
magnetosphere have been calculated assuming incoherent Cerenkov rad- 
iation from magnetospheric electrons with energies from 50 eV to 50 
keV. Calculations were included to determine the ray paths of the 
emitted waves. A diffusive equilibrium model of the magnetosphere 
with an ionosphere, plasmapause, and a centered dipole magnetic field 
was used. Ray path calculations were done in three dimensions. 

Using simultaneous energetic electron and VLF data, comparisons were 
made between calculated and observed VLF hiss. VLF hiss is observed 
to have a sharp lover frequency cutoff near electric and mag- 

netic field spectral densities as large as l.U X lo"® V^/m^z and 
9.6 X lO” Wb /m Hz, and wave normal angles a few degrees (~ 6) from 
the resonance cone. The calculated characteristics of VLF hiss 
include a sheurp lower frequency cutoff near f^ ^ , spectral densities 

of 5 X lO"® V^/m^Hz and 1 x 10”^^ Wb^/m^Hz, and wave normal angles 
“6 “3 

very close (~ 10”° to 10“^ degrees) to the resonance cone angle. 
Assuming a wave normal angle six degrees f^om the resonance cone 
angle, the calculated spectral densities are both two orders of mag- 
nitude below the observed spectral densities. It seems unlikely 
that VLF hiss is produced by incoherent Cerenkov radiation. The 
observed spectral shape of V-shaped VUF hiss is similar to that cal- 
culated from incoherent Cerenkov radiation. However, for V-shaped 
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VLF hiss to have its characteristic spectral shape, the emitting 
region must be small. Calculated spectral densities of incoherent 
Cerenkov radiation from this small region are orders of magnitude 
below observed spectral densities. 

Observed characteristics of ELP hiss have been compared with 
incoherent Cerenkov emd ray path calculations. ELF hiss is> observed 
with spectral densities three and four orders of magnitude higher 
than calculated. Hay path calculations show that Cerenkov. waves ' 
generated at high altitudes may not propagate to the low altitudes 
and over the wide range of latitudes where ELF hiss is observed. 

The observed spectral shape of saucers is similar to that cal- 
culated for incoherent Cerenkov radiation. Calculated spectral den- 
sities are orders of magnitude below observed spectral densities. 

Calculations assinning incoherent Cerenkov radiation of the 
spectral shape and spectral, densities for fast hisslers agree with 
observations. However, »to reach the ground whei;e they have been ob- 
served, they must be generated at initial wave normal angles between 
.0.2° and 2.3° with respect to the magnetic field. Cerenkov ra.dia,tipn 
-is generated .at wave normal angles near the resonance cone. Fast 
hisslers have .very short durations but incoherent Cerenkov radiation 
would be .expected all along the path of. the particles giving longer 
durations . Thus incoherent Cerenkov radiation seems to be i^ed 
out for fast hisslers. 

UHR (upper hybrid resonance) noise is observed to have wider 

spectral shapes and lower spectral densities than caJ.culated from 

.1 
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incoherent Cerenkov radiation. Thermal effects may dominate at these 
ftrequencies and may have to be included in the theory. 

Incoherent Cerenkov and cyclotron emission’ calculations were 
made in Jupiter's magnetosphere. A model of energetic electrons 
emitted and accelerated by lo rather than the electrons in the radia- 
tion belt were used in the calculation. The power levels calculated 
were fourteen orders of magnitude below the power observed in the 
decametric bursts which are associated with the orbital position of 
lo. 



I. IHTRODUCTION 


The magnetosphere is a complex system of energetic particles, 
thermal plasma, and electromagnetic fields. It is clear that the 
sun is the main source of energy for magnetospheric dynamics and 
that the earth and its atmosphere act as the energy sink for the 
system. Understanding the process hy which the energy is transferred 
from the sun to the magnetosphere; how the energy is distributed 
and passed between the magnetospheric energetic particles, themal 
plasma, and electromagnetic fields; and how the energy is lost 
from the magnetosphere is one of the major objectives of the study 
of space physics. In particular, the study of the interactions 
between the energetic particles and the fields in the magnetosphere 
is thought to be the key to the complete understanding of the 
dynamics of the magnetosphere. There are a number of possible wave- 
particle interactions that may be of importance in the magnetosphere, 
and, to determine which dominate, quantitative studies of the 
theories must be done to compare their properties to magnetospheric 
data. It is Just such a study which has been undertaken here to 
determine whether Cerenkov emission is responsible for the existence 
of ELF, VLF, and LF electromagnetic waves in the magnetosphere . 

The Cerenkov mechanism has been a popular suggestion for 
explaining wave observations. Its popularity has been based on 
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some of its in5)ortant characteristics: 

( 1 ) It is an incoherent mechanism so that it can occur regardless 
of plasma stability conditions. 

(2) The emitted power comes primarily from the copius low-energy 
particles. 

(5) The emission is broadband covering the frequency range of 
observed radio noise. 

Previous calculations have had, at best, marginal success in explain- 
ing observations. The calculations to be presented here tend to indi- 
cate that incoherent Cerenkov emission may be responsible for UHR 
noise but is probably not responsible for VtF hiss, fast hisslers, 
and saucers. 

In section II the theory of Cerenkov emission is reviewed, 
and a formula for the power emitted from a particle in a magneto- 
plasma including collisions is derived which, in the case of zero 
collision frequency, agrees with other derivations in the literature 
which were done for the case of no collisions. Section III describes 
ray tracing programs developed from Shawhan [1966, I967] and the 
method of calculating electric and magnetic fields along the ray 
paths. Sections IV through VIII apply the results of sections II and 
III to particle and field data. Discussions of possible problems of 
the theory, including its possible inapplicability in some cases, and 
other effects and mechanisms which might be more successful in ex- 
plaining some wave observations are given in section IX. Section X 
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gives results of Cerenkov calcultions done in the Jovian magneto- 
sphere. Section XI summarizes the results of the Cerenkov and ray 
path calculations. 



II. CERENKOV AND CYCLOTRON RADIATION 


IN A PLASMA 

A. Cerenkov Radiation 

A particle •will produce Cerenkov radiation at a frequency f 
in passing through a medium ■whenever the velocity of the particle 
exceeds the phase velocity v^Cf) of a wave in the medium. To see 
how ,a particle can produce Cerenkov radiation consider the case of 
a test electron moving in a straight line in a plasma. A mo-vlng 
electron will interact with the ambient electrons ajid ions of the 
plasma, primarily causing a momentary imbalance of the local charge 
neutrality by pushing electrons out of the vicinity and attracting 
ions to the vicinity of the moving electron. These interactions 
will therefore cause a momentary charge-separation electric field 
and a momentary current-induced magnetic field which results in a 
small electromagnetic pulse. The frequency components of the pulse 
which have phase velocities faster than the particle interfere 
destructively so that, at a distsint point, the resultant field 
intensity is zero. The frequency components of the pulse which 
travel slower than the particle add coherently to produce a radiation 
field at an angle 9 to the direction of the electron. 9 is given by 
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where n is the index of refraction of the medium at frequency f, 

V is the velocity of the particle, and c is the speed of light. 
Besides emitting Cerenkov radiation, a particle in a magnetic field 
with a component of velocity perpendicular to will radiate 
cyclotron radiation at an angle 9 given by 


u). 


cos 0 = 




ige 


where 

s = ± 1, ±2, ... ; 

u) = 2TTf, f = frequency; 

U) = electron gyrofrequencyj and 

gs 

cPg ~ '^2 ~ velocity parallel to 3^. 

B. Theoretical Treatments of Cerenkov Radiation 
in a Magnetoplasma 

The two types of calculations that have been used to analyze 
the radiation from a particle gyrating in a magnetoplasma are the 
Hamiltonian method and the Fourier transform method. In the 


Hamiltonian approach the vector potential is expanded in a Fourier 
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series, which transforms a complicated set of partial differential 
equations for the vector potential into a system of ordinary differ- 
ential equations for the coefficients of the series. Eidman [1958» 
1959]^ Liemohn [ 1965 ]^ and Trulsen and Fejer [1970] used this method 
to calculate Cerenkov and cyclotron power in the absence of collisions. 

In the Fourier transform method the electric field is trans- 
formed assuming an infinite series of plane waves. This procedure 
allows expansion of two time-dependent factors in a series of Bessel 
functions which reduces the problem to an integral of a series of 
Bessel functions. Mansfield McKenzie [19^7 ], Melrose [I 968 ], 

and H. C. Ko [private communication, 1972] have used this method to 
derive an expression for power in a cold, collisionless magneto- 
plasma. Trulsen [1971] has extended this method to include thermal 
effects through a temperature-dependent dielectric tensor. The 
resulting expressions frcan both approaches are in agreement. The 
derivation outlined here follows Mansfield's but includes the effects 
of collisions and corrects an error in one of Mansfield’s symmetry 
relations. 

The power radiated by an electron in moving through a plasma 
can be calculated by equating it to the work per unit time done on 
the particle by the charge-separation electric field. If q is the 
charge, F is the force, v is the velocity of the electron, and E 
is the electric field at the electron, 

-+ -¥ 

F = qE , 
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P = ? • 


V 

q 


= E • 



( 1 ) 


or, in terms of the current J (= 


P 



Equation (l) will be evaluated to determine the power emitted as 
Cerenkov radiation. 

To determine graphically if Cerenkov radiation is possible for 
a set of plasma parameters and particle energies, the index of re- 
fraction surface may be used. The index of refraction surface is the 
locus of the tip of a vector n = (c/(u)k, where k is the wave normal 

vector. The surface is longitudinally (rotationally) symmetric about 

-¥ 

the static magnetic field and so is u=!ually represented in two 
dimensions. Figure 1 shows a representative index of refraction 
surface and the procedure for using it to determine the existence 
of Cerenkov emission for a particle with a given velocity parallel 
to the static magnetic field. On the index of refraction surface 
for a frequency f a line is drawn perpendicular to at a distance 
l/pg above the origin, where cp^ is the velocity of the particle 
along B^. If the curves intersect, the Cerenkov condition is 
satisfied and the particle will radiate. For Figure 1 it is easiest 
to see the details in the schematic representation. In this particular 
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case the component of n along does have a local maximum at 9=0 
(as shown in the logarithmic representation). Consequently waves 
with four different wave normal angles will be radiated by that part- 
icle. A qualitative characteristic of Cerenkov radiation in the 
magnetosphere is that it is generated with wave normal angles close to 
the resonance cone angle where the index of refraction can exceed the 
longitudinal index of refraction for several orders of magnitude. 

In a collisionless, cold plasma there are no dissipative forces, 
and the index of refraction approaches infinity as 9 approaches the 
angle called the resonance cone angle. That means that a particle 
with a very small energy, even thermal energy, could Cerenkov radiate. 
As will be shown later, for small velocities the energy radiated from 
an electron is proportional to l/v, so to keep the energy radiated 
finite, as it must be, a limiting process is included in the theory. 
Limiting processes which might be evoked are: 

(1) energy limits (a particle may not radiate more energy than it 
has) and 

(2) index of refraction limits (by collisions or thermal effects). 

In the derivation to follow, collisions have been assumed to be the 
limiting process. Collisions are introduced into the formalism 
throxigh a phenomenological collision frequency. 

C. Derivation of Cerenkov and Cyclotron 
Radiation Formulas Including Collisions 
The assumptions made in deriving the Cerenkov and cyclotron 
emission expression are the following: 
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(1) The plasma is cold, homogenous, non-permeable (|j. = and 

immersed in a homogenous and static magnetic field. 

(2) The collision frequencies are much smaller than the wave fre- 
quencies of interest. 

(5) The test charges and the electrcmagnetic waves they emit can 
be neglected in the description of the plasma. 

(4) Collisions are included to limit the index of refraction. 

( 5 ) The energy emitted by the test particle is much less than its 
kinetic energy. 

As Eq. ( 1 ) states, the power emitted frcm an electron is 



The equation of motion including a collisional drag force is 


qj(E + Vj X B^) 


m. 

3 




V 

d 




where q., m., v., and u) . are the charge, mass, velocity, and 
J J J ^0 

collision frequency, respectively, of the jth specie, E is the 
electric field at the particle, and is the static magnetic field. 
Using Maxwell’s equations, a formula may be derived for the Fourier 
transform of the electric field 




( 2 ) 
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Here J is the Fourier transform of the current due to the test 
<1 

charge, and ? is a tensor which includes all of the properties of 
the plasma and The determinant of T set equal to zero is the 

dispersion relation for waves in the plasma. To find the power 
the inverse Fourier transform of the electric field, as given in 
Eq. (2), dotted into the velocity must he evaluated: 


P = q 


2(g, ») 



(5) 


To calculate cu) the test charge velocity is expanded in a 

set of J Bessel functions. The integral required to Fourier 

/“4 \ /“^ \ 

transfoim J (r, t) to J (h, (jo) has a delta function in frequency 
as a factor. 


6(s tUgg + nu) cos e 


where 

s is the order of the Bessel function, 

u) is the gyrofrequency of the test electron, 
ge 

n is the index of refraction, 
to is the wave frequency, 

9 is the wave normal angle = < k), and 

cPg is the velocity of the particle along 
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The integration over frequency results in a non-zero power only when 


J is non zero or when 
<1 


s u) + ncu cos 0 - 0 ) = 0 

g0 ^ 


or 


cos 


0 = 






(4) 


Notice that for s = 0, cos 0 = which is the Cerenkov condition. 

The other s*s are identified with cyclotron emission: s > 0 being 

normal cyclotron emission and s < 0 being anomalous cyclotron 

emission. The time -dependent factor in Eq. (3) is also expanded 

in a series of J Bessel functions, and the orthogonality of Bessel 

functions is used to convert the product of two series of J Bessel 

functions to a series of products of J Bessel functions. The 

integral over ic in terms of J 's may be converted to one over n 

s 

and integrated by contour integration to an expression containing 

products of J's and = J + iY I, where J and Y are 

US s sj s s 

the Bessel functions of the first and second kind of order s. 

At this point in his derivation Mansfield made an error in the 
symmetry relations {Eqs. (55) in Mansfield [196?]} which allowed him 
to incorrectly discard two terms. The correct symmetry relation for 
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the complex calculation including collisions causes the imaginary 
parts to cancel in the frequency integral. Making the correction, 
the final expression for the power spectral density is 


df 


= Re 


r -tt q.^(-i)^ f y j' K' ’V 
•S j ^ 1 2 jr L [\ Jg ”s ^11 

J=1 6. e, c n, - n„| s=-«» 


"o "l ^2 


2 

"2 


+ V. 


T 

s 


H 


( 1 ) 


"22 


+ v" J H 
2 s s 


( 1 ) 


^33 


2v^ -2- 


J H 
s s 


( 1 )' 


"12 


2ViV2 




( 1 )' 


"13 


aviVg 


~ J 
L s 
s 


H 


( 1 ) 


"23 


n=n , 


(5) 


where 

f is the frequency of the emitted wave = ti/2rr, 

Sq is the permittivity constant, 
c is the speed of light, 

Vj^ is the component of the particle velocity perpendicular 
to 

Vg is the component of the particle velocity parallel to 

J is the Bessel function of order s, 
s ' 
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r(l) 


H' is the first Hankel function of order s 
s 


H 


= J + iY 
s 


Y is the Neuman function of order s, 
s ^ 

L is the argument of all of the Bessel functions (J, Y, H), 
s 


00 


L = P- n sin 0 , 

ge 


primes mean the derivative with respect to the argument, 

/-► -4 . 

0 is the wave normal angle [0 = < (k, B ;J defined by 

s so 

Eq. (U), 

s is the mode number, 


2 . 2 


2 2 


T^l = 62.^5 “ ^ 1 " 9<= “ 9 = > 


s "3^ 


24 2 2 

T22 = - e^n + (n - e^n ) sin 0^ 


2 2 2 4 2 2 

T55 = " ^2 " + (n - e^n ) cos 9^ , 


Ti 2 = sin^ 0^ - , 


^15 ^ 2 “^ ®s ®s " 


/ 4 2\ 

= (n - €2_n ) sin 0^ cos 0^ , 


^1 = 1 


-E^ 


2 2 


J CD “ (U, 


fed 
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is the sign of the Jth specie. 




Re means to take the real part of. 


ln=n, means to evaluate the expression inside the bars 
J 

for a particular index of refraction n.. 





22 


C = 
n 


/ “fee 

1 - s 
UL 

V 


'1 ^2 




+ €, 


=1 ^2 


00 . is the phenomenological collision frequency of the jth 
uj 

specie. 


“fee “^ee “^ei * 


‘"ci = “'ii ‘"ij ‘"ie > 

where and cu^^, the s'elf-collision frequencies, are given by 
Spitzer [1956]; ujg^, the electron-ion collision frequency, and 

the ion-ion collision frequency, are given by Nicolet [1953]; 
and the ion-electron collision frequency, is given by Sachs 

[1965] in terms of cu^g. Note that every quantity is complex except 
q, tJD# c, g, v^ g, (jo^^., and s, and that this is due to the 
addition of collisions which adds an imaginary part to the plasma and 
gyrofrequencies. Mansfield's original calculations left out the 
and the ^15 terms for his expression 'for s / 0, but did lead to 
the correct expression for s = 0 and terms). 

Mclienzie [1967], Melrose [I968], and Trulsen and Fejer [-1970] have 
all published results which agree with Eq. (5) for the collisionless 
case. H. C. Ko [private communication, 1972] has also redone 
Mansfield's calculations and agrees with these results. For the 
details of the derivation the reader is referred to ^fa.nsfield [I967]. 
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For s = 0 and small velocities, the power radiated by the 
term dominates. It is proportional to l/v and the others vary as v. 

A program has been developed which calculates the power per 
unit volume emitted by a spectrum of energetic electrons at a fre- 
quency f for specified plasma parameters. Eq. (5) is evaluated for 
each of eighteen pitch angles, 2.5°, 7-5°> •“> 87 - 5 °^ for each 
energy included in the electron energy spectrum. The energy spectrum 
is approximated by specifying the equivalent number density of elec- 
trons in energy intervals centered on a finite number of energies. 

The program calculates the number density of the energetic electrons 
in each pitch angle segment (electrons/dV)„ ^ , and multiplies it by 

Ei,aj 

the power emitted per par*ticle per unit frequency (dP/df electrons) 
to get the power emitted per unit volume per unit frequency 
(d^P/dVdf)„ . The s\im 

i!,i,aj 



is then performed. A representative wave normal angle 0 [corres- 
ponding to the largest (d P/dfdV)„ _ ] is then chosen to represent 

£i* •Cx • 

the direction of the wave normal angle for d'^P/dfdV. (in general, 
the spread in 9 values is 1° . ) Two cards are punched for input 
to the ray tracing program, each with one half of d P/dfdV as the 
volume emissivity and with 9 and -6 as wave normal angles. This means 
that the power is not assumed to be radiated in a cone, but in two 
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directions, corresponding to the intersection of the cone and the 
magnetic meridian plane. 
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III. RAY TRACIKG AND FIELD CALCUIATION 


A. Ray Tracing Routines 

After the power is calculated from a spectrum of electrons it 
is desirable to determine the path the emitted power follows throxigh 
the magnetosphere. Shawhan [1966, 19^71 has developed a computer 
routine to do this calculation in the meridian plane. Using this 
routine as a basis four routines were written, including collisions 
in the ray tracing expressions and including a calculation of the 
coUisional damping the wave experiences as it propagates. Two of 
the routines were for two spatial dimensions (r, 9; i.e., in the 
meridian plane) and two were for three spatial dimensions (r, 0, cp). 
The dependent variables in the two-dimensional routines were, in 
the first case, r, 0, and t where r is the radial distance from the 
center of the earth, 9 is the colatitude, and 'I' is the wave normal 
angle; and, in the second case, r, 0, p^, and pg where p^ and Pg are 
the components of the index of refraction in the r and 0 directions. 
The dependent variables in the first three-dimensional routine were 
r, 0, cp, A, and e where cp is the longitude, A is the angle between the 
radial direction and the vector k (the projection of the k vector in 
the meridian plane), and e is the angle between k^ and k. In the 
second the dependent variables were r, 0, cp, p^, pg, and p^ where p^ 
is the component of the index of refraction in the cp direction. 
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All of the routines were checked by reproducing results of 
Kimura [ 29 ^ 6 ] and by consistency checks. Validity checks will be 
discussed further in section IX. The versions using angles were 
considerably faster computationally because only one (two) angles 
in two (three) dimensions are necessary to completely specify the 
index of refraction. Two (three) components are required in the 
versions using the p's. In the over-determined routines (versions 
using p's) the index of refraction must be corrected at each inte- 
gration step which decreases integration accuracy and hence decreases 
execution speed. Consequently the routines using angles were used 
for all of the results given here. 

B. Electric and Magnetic Field Calculaticns 
To compare the theory to data it is often desirable to calcu- 
late spectral densities of fields from the spectral density of power. 
Mosier and Gurnett [1971] have already made the conversion from fields 
to powers in their Eqs. (l), (2), and (4). They are as follows: 


E = e cos 

(-U*) , 

B = b 

sin 

(-out) 

X X 


X 

X 


E = e sin 

(-ut) , 

B = b 

cos 

(-out) 

y y 

y 

y 


E = e cos 
z z 

(-ut) , 

B = b 
z z 

sin 

(-out) 


where 
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e = E . 

X o ^ 


b = 

X 


n cos 0 2 


s - n 


2 o 


V _ n 


2^0 ^ 


b = 

y 


n cos 9 




6^ - n sin 9 


2 € 

^ _ n cos 9 sin 9 ^ _ ti sin 9 2 „ 

z 2. 2. o' z c 2o 

e^ - n sin 9 €^ - n 


E. = <sj 




z' n cos 9 


2 . 2 ^ 
e^ - n sin 9 


(^1 - n 


2V 


(s ) is the z component of the time averaged Poynting flux, 
z 


<S^> = S cos cp, cp = < (^Q>^g) =0+0: , and 


a = < (i^, ^ ) = tan-^ (- 

' ' g'^ I n de 


0 = < (k, B^) . 


The geomagnetic field is in the z -direction and the wave normal 
angle for the two-dimensional progreim is in the x-z plane. The cal- 
culation of the electric and magnetic field components is done 
point-by-point in the ray path calculation. 
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V 5 IV. . VLF HISS 

' ) ' '1 

A. Observations 

VLF hiss has been described by a number of investigators 
[Helliwell, 19^5 j Gumett, I966; Jorgensen, 1966; McEwen and 
Barrington, 19^7; Laaspere et al., 197I; Gurnet t and Frank, 1972; 
and others]. Laaspere et al. [I971], Barrington et al. [1971], 
and D. A. Gumett [private communication, 1972] chajracterize the 
spectral shape of VLF hiss in low-altitude satellite observations 
as having a sharp lower-frequency cutoff, the largest spectral 
density just above the lower-frequency cutoff, and an intensity 
which is slowly monotonically decreasing with increasing frequency 
above the maximum spectral density. The lower-frequency cutoff is 
generally thought to be the local lower hybrid resonance frequency 
where downgoing waves are reflected so that lower frequency waves 
propagating down are not observed. Jj^rgensen [1968] reported low- 
altitude satellite observations of VLF hiss which seemed to have a 
maximum near 10 kHz with gradually decreasing fluxes above and below 
10 kHz. However, he drew this conclusion from 26 two-point spectra, 
all of \rtiich could also be fit with the spectra found by Laaspere 
et al. [1971], Barrington et al. [I97I], and D. A. Gurnett [private 
communication, 1972]. Kie lower- frequency cutoff of VLF hiss is 
usually 2 to 6 kHz with upper- frequency cutoffs as high as 5^ kHz. 
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Figure 2 shows an Injun 5 VLF hiss spectrogram taken from Hosier and 
Gumett [1972]. 

A low-altitude satellite usually observes VLF hiss for only 
a matter of minutes, but this variation is considered to be spatial 
rather than teniporal since observations of VLF hiss have been made 
on consecutive orbits ninety minutes apart at the same location, 
indicating diirations at least occasionally of a few hours. 

Gximett [1966], using Injiin 5 data, found that the spatial 
occurrence of 5.5— 7*0 kHz VLF hiss ^ 5 X 10~ 7 /Hz was predomi- 

nately from 12 to 24 hours magnetic local time and 7“ wide in invar- 
iant latitude, with the latitude of peak occurrence being about 77° 
at l4 hours decreasing to 70° at 22 hours. 

Gumett and Frank [1972], using Ihjun 5 data, found an associa- 
tion between VIjF hiss in the auroral zone and high fluxes of low- 
energy electrons. They found the associated VLF hiss and electrons on 
field lines connected to the dayside polar cusp and to the downstream 
magnetosheath and distant plasma sheet which are the topological 
extensions of the polar cusp into the magnetotail. Gumett and Frank 
also found the VLF hiss to be generally propagating down the geomag- 
netic field lines with power fluxes occasionally as large as 
10 10~^ w/mSiz at 2.5 kHz. (See also Hosier and Gurnett [1972].) 

The electrons associated with these VLF hiss events had energies on 

Q 

the order of 100 ev to several kev and fluxes of about 10^ electrons 

p ^ *1 

(cm sec ster)~ . In addition, a ’threshold' effect was found 
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between the VLF and electron fluxes. For electron fluxes below 
10^—10^ electrons (cm^ sec ster)”^ at 100 ev little or no VLF hiss 
was observed, but above this ’threshold' Gumett and Frank found that 
VLF hiss was essentially always observed. However, VLF hiss was 
found to be not necessarily proportional to the low-energy electron 
flux; In one case (Figvire 4 [Gurnett and Frank, 1972]) 
dj/dE(l60 s E s 280 ev) increased one and one-half orders of magni- 
tude, but dE /df (f = 7«55 kHz) increased three and one-half orders 
of magnitude. 

Since VLF hiss is predominately downgoing, G\irnett and Frank 
[1972] suggest that the generation region is above the satellite, 
and have made a rough estimate of the altitude of the generation 
region to be about 5000 to 10,000 km. 

Laaspere et al. [1971] have reported an estimate of peak flux 
for VLF hiss of 6 X lO” w/m^Hz at 'audio' frequencies from OGO-6 
data. This agrees irtth the measurements of Gurnett and Frank [1971] 
even though OGO-6 had only an electric field detector. To calculate 
power fluxes from electric fields Laaspere et al. [1971] assumed that 
the propagation was longitudinal, that the wave normal angle was zero, 
and that the index of refraction was unity. 

With Alouette 2 electric-field data Barrington et al. [1971] 
reported an average signal intensity of 2 x 10~^ w/m^z at 10 kHz 
assuming an index of refraction of 10 and an average over a range 
of wave normal angles. 
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B. Mechanisms for Producing VLF Hiss 
Various mechanisms have been suggested for the generation of 
VLF' hiss. The Cerenkov mechanism has been considered by Ellis 
[ 1957 , 1959 , i960], Dowden [i960], McKenzie [1963, 196 ?], Liemohn 
[1965], Seshadri and Tuan [1965], Gershman and Trakhtengerts [1966], 
J/iJrgensen [1968], Hartz [ 1970 ], Trulse'n and Fejer [ 1970 ], Idm and 
Laaspere [ 1972 ], and James [ 1973 ]. Another popular mechanism sug- 
gested for VLF hiss is cyclotron emission, described by Ellis [ 19591 , 
MacArthur [ 1959 ], Murcray and Pope [i960, I96I], Santirocco [i960], 
Gintsburg [1961], Gershman and Trakhtengerts [1966], Trulsen and 
Fejer [ 1970 ], and Trulsen [ 1971 ]. Another mechanism suggested is 
a plasma instability, the traveling wave tube mechanism of Gallet 
and Helliwell [ 1959 ] and Dowden [1962]. 

C. The Cerenkov Calculations of Jj^rgensen 
and Lim and Laaspere 

Of the Cerenkov calculations, Jji^rgcnsen ' s [1968] and Lim and 
Laaspere ’s [ 1972 ] are the easiest to compare to VLF hiss data. 
Jjf^rgensen used Mansfield's [I967] theoretical results for the power 
emitted by a particle gyrating in a magnetoplasma and calculated a 
VLF hiss spectrum he claimed would be observable in the lower iono- 
sphere. Jjirgensen used an intense but reasonable electron spectrum 
determined by Evans [1966] for energies above 1 kev shown in Figure 5 
and he assumed an isotropic pitch angle distribution. He also assumed 
the power produced to be independent of the pitch angle of the 
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energetic electrons. He chose the field line corresponding to 70° 

invariant latitude (L = 8.5) and assumed perfect guiding of the 

emitted waves. The guiding allowed him to add all the power produced 

in a flux tube. He chose the upper limit for power production to be 

26,000 tan. Jorgensen Justified his assun^Jtion of perfect guiding on 

the basis of ducting. . J}(5rgensen 's results are shown in Figure 4. He 

-l4 / 2 

calculated a maximum spectral density of about 10 w/m Hz at about 
10 kHz. This is two to three orders of magnitude less than the most 
intense VLF hiss spectral densities detennined by Gurnett and Frank 
[1972] and Hosier and Gurnett [1972] using Injun 5 data, and two 
orders of magnitude below the maximum VLF intensity measiired with 
Injun 3 [Gurnett, I966]. In the same paper Jjf^rgensen reported OGO-2 
observations of VLF hiss. The maximum spectral density was 
4.4 X 10" w/m Hz at l8.5 kHz with a typical average spectral density 
of 10~^^ — 10~^^ w/m^Hz. The spectral densities reported were at 
frequencies from 2.3 to 88 kHz. 

However, j/rgensen's [I968] OGO-2 power spectral densities 
were obtained from the observed magnetic spectral densities using 
a formula from Helliwell [I965] which is applicable only for longi- 
tudinal propagation (wave normal angle = 0° ) . As shown qualitatively 
in section II, for the Cerenkov and cyclotron mechanisms the wave 
normal angle is considerably different from zero and is near the 
resonance cone angle where the index of refraction n is large. This 
large n allows the phase velocity of the wave to be comparable to 
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the velocity of the radiating particle. This condition is necessary 
for most wave-particle interactions. It is likely that the emitted 
wave normal angle of the VLF hiss wave is close to the resonance cone 
no matter what the mechanism. Ray path calculations of waves emitted 
above satellite altitudes with wave normal angles near the resonance 
cone show that wave normal angles stay near the resonance cone. 

These paths are shown in Figure 5* Therefore it seems that 
J;^rgensen's conversion from magnetic fields to powers is not justi- 
fied. If Jj^rgensen ' s observations were of waves with non-zero wave 
normal angles, the power calculated was too low since magnetic fields 
get smaller as the wave normal angle gets close to the resonance 
cone angle. 

Also, a second assumption that Jorgensen made does not seem 
to be justified. He assumed that the waves were perfectly guided 
from 26,000 km to satellite altitudes. However, ray tracings of 
Cerenkov generated rays indicate that rays with frequencies above 
2 kHz are dispersed over a larger latitudinal and longitudinal extent 
than the particle precipitation region. Consequently, this dispersion 
leads to a power loss factor up to a factor of 100. Representative 
dispersed ray paths for 2 and 10 kHz from 6250 and 8750 km are shown 
in Figure 6. 

The most recent calculations of VIF hiss, by Lim and Laaspere 
[1972], have made improvements over the calculations by Jj^rgensen 
but still fall short of ejfplaining the maximum VLF hiss intensities 
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observed. They also used Mansfield^ s formula for Cerenkov power. 

They made the proper corrections and recognized the problem of infi- 
nite power radiated from low-energy particles at resonances, and 
limited the index of refraction to 400. The electron spectriom chosen 
by Lim and Laaspere was similar to Westerlund’s [1968] but softer. 
(See Figure 5. ) Isotropy was also assumed. Calculations were made 
at frequencies from 4 to UBO kHz and, as in J^^rgensen [I968], they 
added all the power in a flux tube. They chose an invariant latitude 
of 78°. Their calculated hiss spectrum is also shown in Figure 4. 

Two orders of magnitude still separate the maximum observed Injiin 5 
and Lim and Laaspere 's calculated VLF fluxes. 

To avoid the power infinity Lim and Laaspere placed a limit of 
400 on the index of refraction, corresponding to the condition that 
the phase velocity of the emitted wave was twice the thermal velocity 
of the background electrons for their assumed electron tenperature of 
5000 °K. However, Stix [1962] and Kennel and Petschek [1966] give 
the analog of Laudau damping for oblique waves to be when knv|| = cw. 
Physically the guiding centers of the plasma particles are confined 
to move along the magnetic field, so that the interaction is between 
the projection of the phase fronts onto the direction of the geomag- 
netic field, and the resonant particles which are constrained 
to first order motion along For example, at aji altitude of 

10,000 km and a geomagnetic latitude of 70° the resonance cone angle 
at 4 kHz is about 80°. An n || of 400 corresponds to an n of about 
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2200. In other words, Lim and laaspere should have placed the limit 
on nil rather than on n. Restricting n|| could make an order-of- 
magnitude increase in their power calculations which would be in 
agreement with the upper limit calculation presented in the next part. 

Limiting n rather than ny probably increased Lim and Laaspere 's 
frequency of peak power. Slightly above the resonance cone 

angle is close to 90°^ and although n may be very large, n|| is small, 
and the power at these frequencies was not included by Lim and 
Laaspere. Since the region of resonance, and thus high Cerenkov 
power, in the frequency range is < f < min (fgg> ^pe^^ 

Laaspere -underestimated the power at the lower frequencies. 

The spectral shapes of VLF hiss calculated by j/rgensen [1968] 
and Lim and Laaspere [1972] were dissimilar. J/^rgensen's was peaked 
at 10 kHz, with a gently decreasing spectral density above and below 
10 kHz. This spectrum -was consistent with the VLF hiss spectrum 
he deduced from twenty- six obser-vations of spectral density at two 
frequencies from the set [2.3^ 5.1, 8.4, 11, l8.5^ ^0^ 60, 88 kHz]. 

It is not consistent with the VLF hiss spectra observed by Laaspere 
et al. [1971], Barrington [1971]^ and D. A. Gurnett [pri-vate communi- 
cation, 1972]. Lim and Laaspere 's spectrum had a maximum at 70 
kHz, gently decreasing below 70 kHz and dropping off rapidly above 
70 kHz. This spectrum is considerably different from any VLF hiss 
spectrum thought to be typical. ' 
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D. Cerenkov Production of VLF Hiss 
1. Upper Limit Calculation 

Low-energy electrons which have large fluxes in the auroral 
zone will emit significant amounts of energy in the Cerenkov mode 
only for frequencies which may propagate with a large index of refrac- 
tion. The frequencies for which this condition holds in the VLF hiss 
range is f^ < f < min (f^^, f^^). It is in this frequency range 
that calculations of incoherent Cerenkov noise similar in spirit to 
Jorgensen's [1968] and Lim and Laaspere's [1972] have been made. An 
even more intense electron spectrum measured with an electrostatic 
analyzer on the Twins I rocket in an aurora^ from Wester lund [1968], 
was used. (See Figure 5. ) This spectrum is somewhat more intense 
than the ones measured by Injun 5 LEFEDEA's [Ackerson and Frank, 1972] 
during an ~ 3 X 10~ w/m Hz VLF hiss event described by Mosier and 
Gurnett [1972]. The maximum electron flux is 1 to 2 orders of magni- 
tude larger than the one used by jOrgensen and the spectrum is ex- 
tended down to about 45 ev. Westerlijnd's spectrm was the most in- 
tense that could be found in the literature. Isotropy was assumed, 
and power was calculated in each of eighteen pitch angle segments. 

The index of refraction was limited by including collisions and by 
requiring that the phase velocity of the emitted wave along the 
magnetic field direction be higher than the thermal velocity of the 
background electrons. For an electron ten5)erature of 5000 °K this 
limits nil to about 6OO. This procedure avoids the problem of infinite 
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calculated power. The background plasma was represented by an 

electron-three ion diffusive equilibrium model with an ionosphere 

having a maximum electron density of 10^ cm ^ at 150 km altitude and 

a plasmapause at L = 4. Two altitude profiles from the model are 

shown in Figure 7 for magnetic latitudes of 0° and 70°. A centered 

dipole was used as a model for the earth's magnetic field. The 

calculations were made at 70° invariant latitude and all the power in 
2 

aim flux tube was added assuming perfect guiding. Upper-altitude 
limits were chosen to be 26,000 km (as Jj^rgensen [I968]) and 48,000 
km (the equator). Account was taken of the flux tube e:q)ansion. The 
results are shown in Figures 4 and 8. Figure 8 shows power produced 
in 2000 km altitude segments for different frequencies and Figure 4 
shows the power summed along the field line. The power peak of about 
7 ^ 10~^^ w/m^z at 10 kHz is nearly as large as the power estimated 
at 8.8 kHz from In,1\m 3 data by Gumett [1986]. The calcrilated power 
at 5 kHz, 2 — 5.5 ^ 10”^^ w/m^z, is two orders of magnitude below 
the power of about 5 X 10“^^ w/m^z measured by Mosier and Gumett 
[1972]. As with previous calculations, the calculated spectrum has 
a maximum at the wrong frequency and has an incorrect lower- frequency 
cutoff. 

This type of calculation gives only a rough estimate of the 
maximum incoherent Cerenkov power which might be observed at satellite 
altitudes. A more accturate calculation would determine the exact ray 
paths of the waves after they are emitted by the Cerenkov mechanism. 
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The waves are approximately guided by the earth's magnetic field but 
do diverge in latitude and longitude. The energetic electron spectra 
used should be that measured at the altitude of generation and it is 
not ejqjected to be constant as a function of altitude. The attenua- 
tion of the waves as they propagate from the generation region to the 
observation region should be included. Experimentally, wave power is 
calculated from one or more field coitiponents, whereas in a theoreti- 
cal calculation the fields can be determined exactly from the power, 
wave normal angle, and plasma parameters. It seems reasonable then 
to compare the directly measured quantities, fields, with the cal- 
culated fields rather than a derived quantity, power, with calculated 
power. 


2. Detailed Calculation 

In order to make a more accurate estimate of what VLF hiss 
intensity a low-altitude satellite would observe from electrons 
radiating in the Cerenkov mode, a more sophisticated calculation has 
been performed. Cerenkov power was calculated from an electron spec- 
trum measured during the same VLF hiss event that was characterized 
as among the most intense observed with Injun 5, shown in Figure 2 . 
Mosier and Gurnett [1972] reported VLF hiss observations and Ackerson 
and Frank [1972] reported plasma measvirements during the same auroral- 
VLF hiss-electron precipitation event on December 21, I 968 . The 
energetic electron spectrum reported was used to calculate Cerenkov 
power at altitudes up to 10,000 km, which Gxxrnett and Frank [1972] 



suggest might be the maximum altitude of primary auroral electron • 
acceleration. The three-dimensional ray tracing program calculated 
the ray paths from the generation region to an altitude of 2500 km, 
the approximate satellite altitude. 

The locus of points where the ray paths crossed 2500 km was 
approximately circular. The width of the precipitation region was 
estimated from Plate 2 of Acker son and Frank to be about 1.8° of 
arc at the 2500 km altitude level. The power loss factor due to 
coUisional damping and to spreading of the wave energy to a region 
greater than that of precipitation can be estimated to be [(width of 
wave illumination region)/ (width of precipitation region)] X 
[collisional loss]. The power expected at satellite altitude is then 



p 

where d P/dfdV is the power spectral density per unit volume emitted 
by the energetic electrons radiating in the Cerenkov mode and V is 
the volume of the emitting region, the voliime of the ith section of 
the flux tube with one m cross sectional area. The summation was 
done in sections from 2500 km to 10,000 km, with d^P/dfdV, V, and 
loss calculated for each section. Similar calculations were made for 
dE^/df and dB^df. 

Figure 9 shows the spectral densities of the power and of the 
electric and magnetic fields that would be expected dioring the 
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December 21, 1968 event due to Cerenkov radiation. Note that the 
power and fields are very much smaller for frequencies below the 
local Hosier and G\irnett report that the maximum intensity of 

the \TjF hiss occurred at Ih b2tn 4?s and the broadband (0.5 to 10 kHz) 
electric and magnetic field strengths were 8.4 X 10 V/m and 22 mj, 
respectively. From these facts and the facts that the cosine corre- 
lation (cos cp) was measured to be about 0.5 aJid the bandwidth was 
estimted to be 5 kHz, they found the lower limit of the VIjF hiss 
power flux to be 1.5 ^ lo" w/m Hz. Converting the measured fields 
to spectral densities assuming a 5 kHz bandwidth gives 

^ = 1.4 X 10"® V^/m^Hz 


dB r. c ^r ^-26 2/ 4„ 

= 9.6 X 10 Wb /m Hz 

df 


Comparing these spectral densities to Figtire 9 it is seen that the 
calculated spectral density of the electric field is the same order ■ 
of magnitude as the measured one but a little larger, and the calcu- 
lated spectral densities of the power and magnetic field are consider- 
ably below the measured fields by factors of lO"^’^ and lO”^, 
respectively. A combination of two possible effects may accoimt for 
this discrepancy. 
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The first effect involves the dependence of electric and mag- 
netic field amplitudes upon wave normal angle. For a constant 
power, as the wave normal angle of a wave approaches the resonance 
cone angle, the electric field of the wave increases and the magnetic 
field decreases. As mentioned above, the wave normal angle of a 
wave generated by the Cerenkov mechanism is generally very close to 
the resonance cone angle at generation and during propagation. The 
dependence of electric and magnetic field spectral densities upon 
wave normal angle is shown in Figure 10. Figure 10 was calculated 
for the 2500 km altitude level of a 2 kHz wave generated by the 
Cerenkov mechanism at 6250 km altitude. There is a difference of a 

factor of 10^ in the spectral densities between waves with wave nor- 

-6 

mal angles characteristic of Cerenkov production (~ 10 degrees 
from the resonance cone) and waves with wave normal angles only one 
degree from the resonance cone angle. There is some dispersion in 
initial wave normal angles due to the fact that particles with dif- 
ferent energies and pitch angles radiate at different wave normal 
angles. For the situation here, however, this dispersion is usually 
less than 0.1°. If, for some reason (perhaps irregularities in the 
backgroimd plasma or by thermal dispersion), some of the waves pro- 
duced by the Cerenkov mechanism were to deviate by a degree or more 
from the resonance cone angle, the electric field spectral density 
would remain high (from the waves with undeviated wave normal angles) 

and the magnetic field spectral density would increase by a factor 

3 

of about 10 . 
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The in 5 )roved calculation above gives an electric field spectral 
density about 5.5 times the observed one. Therefore only of the 
waves need to have wave normal angles comparable to the ones calcu- 
lated by the ray tracing program. If the other 72^ have wave normal 
angles about a degree from the resonance cone angle, the magnetic 
field measured would be about 720 times the one predicted by the ray 
tracing program, or 

^ = 7.2 X 10"^® Wb^/m^Hz 
df ' 


This spectral density is about two orders of magnitude below the 

observed magnetic field spectral density. Hosier and Gurnett [1971] 
have shown, however, that in the power calculation they perform the 
cosine correlation would be small in a case of wave normal angle 
dispersion for uncorrelated waves. Since the cosine correlation is 
typically near 0.5 for VLF hiss, there must be little wave normal 
angle dispersion. 

The second effect which might account for a higher measured 
than calculated power density is the possible sensitivity of a mag- 
netic loop to electric fields. D. A. Giimett [personal communica- 
tion, 1973] has made a rotigh calculation of the voltage ratio ex- 
pected between the electric and magnetic components of a wave pro- 
pagating in free space for an unshielded 1 meter loop at a frequency 
of 5 kHz. The ratio was determined to be about lO”^. Since the 
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In;jun 5 magnetic loop has 6 turns and has a diameter of O .56 m, the 

relative sensitivity in free space from this calculation is perhaps 

-6 

an order of magnitude less or about 10 . However, the sensitivity 

would be greater for waves in a plasma near the resonance cone since 
the index of refraction would be larger as would the electric field. 
Using the same sample ray calculation as was used in Figure 10, the 

ratio of E to B was calculated and found to be 5-75 X 10^^ m/s, a 

3 

factor of 2 X 10 larger than c, the ratio of E to B in free space. 

In Gurnett's calculation free space propagation was assumed and 
e/b = c was used. A correction for the plasma then is to use 
E/B = c X 2 X 10^ which changes the voltage ratio to ~ 5 X lO”^, 
still an insignificant correction. 

An experimental determination of the sensitivity of the Injun 
5 magnetic loop to electric fields was performed and it was deter- 
mined that at5hHz a 150 mV/m electric field induced a signal in the 
magnetic loop equal to that produced by a magnetic field of 0.1 my, 
the approximate noise level in the loop. The intense VLF hiss event 
described above [Mosier and Gumett, 1972] had an electric field of 
only 8.4 mV/m. It is unlikely, on the basis of this experiment, that 
magnetic field observations from Injun 5 are incorrect. It is con- 
cluded that the apparent magnetic field induced by the observed elec- 
tric fields is below the magnetic receiver noise level and could not 
contribute to the observed magnetic field. 

At the time of maximum measured VLF intensity the e'lectric 
field spectrum of the VLF hiss in Figure 2 has a sharp lower cutoff 




at about 1 kHz, which is similar to the spectioun in Figirre 9- This 
lower cutoff is thought to be at the local so a small change 

in the model would reduce the calculated lower- frequency cutoff to 
agree with the data. The upper- frequency cutoff is not as sharp, 
with noise evident up to perhaps 7 kHz* If the dynamic range of the 
spectral analysis system is about 5 db, this agrees with the cal- 
culated electric field spectrum in Figure 9 which is down to l/2 
its maximum value at 7 kHz. This decrease in power for higher 
frequencies is due to the increased angular dispersion of the ray 
paths as the wave frequency increases. Even thou^ more power is 
generated at higher frequencies the spreading of the power over 
larger areas more than compensates for the increased power, resulting 
in an observed spectrum peaked near the lower cutoff frequency. The 
magnetic field spectrum in Figure 2 appears flatter but still has 
upper- and lower- frequency cutoffs. The calculated magnetic field 
spectrum is also flatter, but it extends to at least 10 kHz. 

Laaspere et al. [1971] estimate the 'audio' frequency power 

-12 

from an intense auroral hiss event to be of the order of 6 X 10 
w/m Hz, which corresponds to an electric field spectral density of 
about lO”^ V^/m^Hz about 50 times lower than the maximum in Figure 9* 
Barrington et al. [1971] estimate an average signal intensity at 
10 kHz of 2 X 10 w/m^z, which corresponds to an electric field 
spectra density of about 5 X 10“^ V^/m^Hz. Jorgensen [1968] esti- 
mated a maximum power spectral density of 4.4 X lO”^ w/m^z 
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corresponding to a magnetic field spectral density of about 10 

2 / 4 

Wb /m Hz, very close to the Injun 5 measurements of Mosier and 
Gurnett [1972]. These calculations and measurements of VtF hiss 
power and field densities are summarized in Table 1. 

The E/B ratio calculated for the event shown in Figure 2 cor- 
responds in the model to an index of refraction of about 8 and a wave 
normal angle about 6° from the resonance cone. Since the cosine 
correlation was foimd to be 0.5, the wave normal angle dispersion was 
at most a few degrees. If the wave normal angle of the Cerenkov waves 
as calculated here were actually 6° fr-om the resonance cone, the 
electric field would be reduced and the magnetic field would be in- 
creased by 3 orders of magnitude. The calculated spectral densities 
in this case would be factors of about 550, 100, and 300 below the 
observed spectral densities for the electric field, the magnetic 
field, and the power, respectively. Even though the spectrums have 
the right shape, the fact that the calculated and observed spectral 
densities are so much different for VLF hiss seems to rule out inco- 
herent Cerenkov radiation as the source of VLF hiss. Other possible 
source mechanisms are discussed in Section IX. 

3. Cyclotron Modes 

The contribution to VLF hiss by other incoherent modes (s > 0, 
cyclotron radiation and its harmonics, and s < 0, anomalous cyclotron 
radiation) was judged to be minor, at least for the most intense 
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lower frequencies. Figure 11 shows the power emitted from the elec- 
tron spectrum used in the above calculation of VLF hiss at an alti- 
tude of 10,000 km and an invariant latitude of 70° for various modes. 
Except very near the plasma frequency, the cyclotron contribution is 
an order of magnitude below the Cerenkov power. 

E. V-Shaped VliF Hiss 

One form of VLF hiss that is observed in low-altitude satel- 
lites is V-shaped hiss, observed to be propagating down the field 
lines, in which the minimum frequency of the hiss first decreases 
and then increases in time, with a total width of 30 seconds or more 
and a minimum frequency of 1 kHz or more.( (See Giomett [1966], 
Gurnett et al. [1971]^ and Gumett and Frank [1972].). Gurnett and 
Frank [1972] have hypothesized that the V-shape is due to the fre- 
quency dependence of the limiting ray angle for whistler mode propa- 
gation from the generation region down to the satellite. In other 

~¥ 

words the maximum angle with respect to at which the wave may 
propagate increases as frequency increases, allowing high frequency 
waves to be observed further from the field line passing through the 
generation region. 

To check the theory of Gurnett and Frank [1972] the V-shaped 

VLF hiss in their Plate 5 (reproduced here as Figure 12) was exam- 

h m 

ined. The first V-shaped VLF hiss event at 22 55 UT in Plate 5 

has a width of about 1.6° in invariant latitude between points where 
the frequency of the outer edge of the V-shape is 5 kHz. The mag- 
netic local time is changing very slowly at this point so that the 
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satellite is moving primarily in invariant latitude. The power pro- 

! 

^ gram was used to calculate initial wave normal angles for Cerenkov 
generated waves for frequencies from 2 — 12 kHz at an altitude of 
4000 km and an invariant latitude of 73°. The ray paths were then 
determined by the two-dimensional ray tracing program. The lati- 
tudes where the ray paths passed an altitude of 2500 km are shown in 
Figure 13 . Although the shape of the expected spectrum from Figure 
13 does not precisely match the shape of the event in Figiire 12, 
their general appearance is similar. The width at 5 kHz is 1.45°^ 
thus indicating that the generation altitude of this event would be 
slightly more than 4000 km if this model is correct. Gumett and 
Frank [1972] estimated the altitude of the VLF hiss generation region 
to be about 5000 — 10,000 km. This is in reasonable agreement, espe- 
cially since this example appears to be the narrowest of the V-shaped 
VLF hiss events in Figure 12. 

The fact that there is less noise inside the outline of these 
events suggests that the characteristic size of the generation region 
is considerably less than the distance between the generation region 
and the satellite. An incoherent Cerenkov mechanism is probably not 
responsible for this VIF hiss since the precipitating electrons 
would radiate over much wider altitude range, producing a more inco- 
herent spectrum. This decrease in power for higher frequencies is 
due to the increased angular dispersion of the ray. Even if the 
particles were confined to radiate in a source region of ~ 100 km. 
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the power radiated would not be sufficient to explain the observed 
power or electric field spectral density. The generation mechanism 
may be an instability or coherent Cerenkov radiation. The wave nor- 
mal angle of the emission is probably near the resonance cone angle 
since the initial wave normal angles for Cerenkov emission were used 
in this calculation and are near the resonance cone. 

If the satellite passed throi:igh the generation field line the 
observed minimum frequency of the emission would be the local 
since these waves would travel with least deviation down the field 
line and reflect (reverse direction) at the satellite altitude. 

Higher frequency waves would be able to propagate below the satellite 
altitude, but lower frequency waves would reflect above the satellite. 
If the satellite were to pass to one side of the field 3jjie through 
the generation region the minimum observed frequency would not be 
the local f that frequency with a limiting ray angle large 
enough to propagate to the satellite. Such V-shaped VLF hiss events 
with variations in the lower cutoff frequency are observed. (See, 
for example. Figure 7 of Gurnett et al. [l97l]*) 



V. SAUCERS 


The last event at 22h 58m UT in Figure. 12 was observed to be 
traveling up the field line. It is an example of a saucer-shaped 
emission, or a saucer. Observed with the Alouette 1 and 2 satellites 
(R. E. Barrington, personal communication, quoted by Mosier and 
Gumett [1969]) and Injun 5 by Mosier and Gumett [1969]^ Gurnett 
et al. [1971], and Gurnett and Frank [1972], saucers are often ob- 
served at low altitudes in the auroral zones. They range in fre- 
quency from less than 1 kHz to above 10 kHz and have a time duration 
from a few seconds to a few tens of seconds. 

Mosier and Gumett [1969] have presented a detailed theory 
explaining many of the characteristics of saucers. Its appearance on 
the spectrogram depends upon the same limiting ray-angle argument 
that was mentioned above for V-shaped VLF hiss. Basically, a source 
region below the satellite emits waves at frequencies above the fj^jp 
at the source which propagate upward. As the satellite passes over 
the source region, high frequency waves with large limiting ray angles 
are observed, followed by lower frequency waves with smaller limiting 
ray angles until the minimum frequency, the s-t the sovirce, is 

observed when the satellite passes through the field line connected 
to the source region. Then the frequency increases as the satellite 
moves away from the source field line and the limiting ray angle 


increases. 
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Figrire lU shows another example of a saucer observed with 
InjTon 5* Mosier and Gumett C1969] have made a detailed analysis of 
this saucer and from its width have deduced from limiting ray-angle 
considerations that the source was located at an altitude of lUOO km, 
about 1100 km below the satellite. Gurnett and Frank [1972] consider 
it likely that intense low-energy (~ 100 ev) streams of electrons 
are the particles which emit saucer radiation. 

It is possible that the generation mechanism of saucers is 
incoherent Cerenkov radiation. To check this possibility and the 
theories of Mosier and Gixrnett [1969] s^nd Gxirnett and Frank [1972], 
Cerenkov power and ray path calculations were performed to match the 
characteristics of the saucer observed at an altitude of 2500 km 
shown in Figure l4. Best agreement between the calculated and ob- 
served characteristics of the saucer in Figure l4, width 12 sec 
at 10 kHz) and lower- frequency cutoff 1.7 kHz), was found for a 
source altitude of iBOO km. Since no energetic electron spectrum has 
been measured in association with a saucer the electron spectrum 
measured by Westerlund [1969] was used in the Cerenkov power calcula- 
tion (assumed to be flowing up the field lines, cc > 90°)* The ray 
paths of the Cerenkov waves were then determined and the latitude of 
the waves at the altitude of 2500 km as a function of frequency is 
shown in Figure 15. The shape of the calculated saucer has no ctirva- 
ture as the observed saucers do. This curvature may be due to the 
small change of plasma parameters in the model over the 700 km ray path. 



53 


Ass\ming an emission region with a characteristic dimension of 
100 km the spectral densities at 8 kHz were calculated to be 


^ = 8 X lO”^”^ w/m^z 
^ = 5 X 10"^^ V^/m^Hz 
^ = 3 X 10"^^ Wb^m^Hz 


The maximum wide band electric and magnetic fields measured during 
the saucer in Figure lU were 5*6 x lO”^ volts and 1.9 ni7- Assuming 
a bandwidth of 2 kHz and an effective dipole antenna length of 2.85 ni, 
the spectral densities can be calculated to be 

^ = 1.9 X 10"^ V^/m^z 


^ = 1.8 X 10"^”^ Wb^/m^Hz 


The satellite electric antenna spectrum analyzer at 7*35 kHz measured 
a maximum voltage spectral density of 1.7 X lO”^ V^/Hz giving an 
electric field spectral density of 
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^ = 2.1 X 10"^ V^/m^z 


which agrees quite well with the measurement derived from the wide 
band electric field. Using the wide band measurements and assuming 
a cosine correlation of (cos cp) = 0.5 (see Gumett and Frank [1972]) 
the lower limit of the power spectral density is calculated to be 


II = 7.4 X 10"^^ w/m^Hz 


The observed spectral densities are considerably higher than the ones 

calculated on the basis of incoherent Cerenkov radiation, i.e., a 
l|. 

factor of 10 in power. Coherent Cerenkov radiation where groups of 
electrons radiate in phase might produce enough power to explain the 
intensity of saucers, or they may be produced by another collective 
action of radiating particles, for example an electrostatic or elec- 
tromagnetic instability. 

Since the envelope of saucers is usually sharply defined the 
emitting region must be a confined region rather than a line, sheet 
or extended source since any of these regions would tend to fill in 
the envelope of the saucer. 

From the results obtained here, saucers appear most likely to 
be generated in a limited region below satellite altitudes at an 
altitude of a few thousand kilometers by the collective action of 
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low-energy upgoing electrons. The radiation apparently must be 
generated at wave normal angles near the resonance cone angle so 
that these ray angles will be large to enable, the waves to diverge 
a few degrees in hundreds of kilometers of path length. 
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VI. FAST HISSLERS 

Siren [1972] reported the observation from the ground VLF 
station at Byrd Station (invariant latitude = 71°) of bursts of 
auroral hiss. Showing evidence of whistler mode dispersion, they were 
of very short duration, a few tenths of a second or less at all 
frequencies. They were observed from 2 to 12 kHz with a maximum sig- 
nal intensity of 

= 1.1 X 10'^'^ w/m^Hz 
df ' 

at 3 kHz. Whistler mode analysis of the fast hisslers by Siren [1972] 
indicated generation altitudes in the range of 10,000 to 22,000 km. 
Siren suggests that Cerenkov radiation at these altitudes may be 
responsible for fast hisslers. 

Waves that are generated by the Cerenkov mechanism from elec- 
trons xri.th energies less than a few tens of kev have wave normal 
angles that are very close to the resonance cone. As these waves 
propagate down the field lines their wave normal angles change, but 
remain within a few hvindredths of a degree of the resonance cone 
angle. When waves near the resonance cone reach the altitude where 
their frequency is equal to the local they reflect and begin 

going up the field line. So if fast hisslers are generated by the 
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Cerenkov mechanism the waves would need to somehow couple to waves 
with smaller wave normal angles to pass the altitude level where 
f = Hissler waves must also make the transition from the iono- 

sphere where the index of refraction, n, is very large to the atmo- 
sphere where n is one. If the transition is thought of as a step in 
n, only waves with values of A, the angle between k and the radius 
vector, close to l80° will be able to refract into the atmosphere. 

A calculation of Cerenkov power was performed for an altitude 
of 15,000 km, an invariant latitude of 71°, and a frequency of 5 kHz 
and the ray paths were calculated for the resulting waves. The 
power observed at a ground station may be calculated with the fol- 
lowing formula: 


- 

dP 

JT 

(a lat) ^ 

wave 

df 

gs 

dfdV 

V 

(A lat) 


T T 
MS lA 


where dP/dfdV is the volume emissivity, V is the emitting volxime, 

(a lat) is the spreading of the wave in latitude (longitude), 

(a lat) ... is the width of the particle precipitation region, 
emitters 

T.,_ is the transmission coefficient due to collisional attenuation 
MS 

along the ray path in the magnetosphere, and is the transmission 
coefficient of the wave resulting from the inefficient coupling from 
the ionosphere to the atmosphere. The energetic electron spectrum 
used was that measured by Ackerson and Frank [1972] during the 

auroral- VLF hiss-electron precipitation event described in section IV. 

I 

I 

I 

I 

i 
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(a lat) was determined with the ray path calculation to be ' 6 . 9 ° 

' wave ^ 

and (a lat) ... assumed to be l.S" as before. The ray path trans- 
emitters 

mission coefficient Tj^g was found to be O .9 (“ 0*^ db). was esti- 

mated from maximum satellite and gro\and station VLF hiss observations 
of 1.5 X 10 ”^^ w/m^Hz [Mosier and Gumett, 1972] and 10 w/m%z 
[Helliwell, I 965 ] to be about lO"^. Knowing the emissivity 
1.07 X 10 w/m^Hz and the observed spectral density 1.1 X 10 
an estimate of 5000 tan for the source region size results. 

One of the fast hisslers presented by Siren [1972] had an 
indicated generation altitude of 15,000 ± 1300 tan. So starting at an 
initial altitude of 15,000 tan it was determined that there was a 
range of initial wave normal angles for which the waves would success- 
fully reach ground level at \ = 71°. These ranges, as a function of 
frequency from 2 — 12 kHz, are shown in Figure I 6 . Therefore, if the 
generation point of this hissler was 15,000 km and the propagation 
of the waves to the ground involved no coupling to other wave normal 
angles, the initial wave normal angles of the waves had to be bet- 
ween 0.2° and 2.5° with respect to the magnetic field. To check to 
see if the model produced the same dispersion as the hissler, ray 
paths for frequencies from 2 — 12 kHz were calculated from an initial 
altitude of 15,000 to a final altitude of 100 km and dispersions com- 
pared to those scaled from the spectrogram of the hissler given in 
Siren [1972]. The comparison is shown in Figure 17 v?ith the line 
drawn through the dispersion determined from the ray path calculations. 
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The s m all difference in the two dispersions indicates that the total 
electron content below 15,000 km was similar in the models used by 
Siren and in the ray path calculation. 

From power considerations, fast hisslers may be generated by 
the Cerenkov mechanism, but if they are the initial waves must couple 
to waves with small wave normal angles to pass the point where 
f = iiiust have k vectors antiparallel to the radius vector 

in the ionosphere to reach the ground. If fast hisslers propagate 
to the ground directly their initial wave normal angles must be close 
to zero. No matter what the generation mechanism the generation 
must be limited to a few thousand kilometers in size and must turn on 
and off on the order of a tenth of a second since the ground observa- 
tions show whistler mode dispersion. 
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VII. UHR NOISE 

UHR noise is band- limited white noise that is commonly ob- 
served by rockets and satellites in the magnetosphere at frequencies 
for which max (f^^, f^^) ^ f ^ where f^ = f|^ + f^^. 

Reports of upper hybrid resonance (UHR) noise in the magneto- 
sphere have been given by Walsh et al. [1964] using rocket data, Bauer 
and Stone [1968] using ATS III data, Harvey [I968] using Ariel II 
data, Gregory [1969] using Ariel III data, Barrington and Hartz 
[1969] using Alouette I and II data, Hartz [1969> 1970] and 
Barrington et al. [1971] using Alouette II data, Muldrew [1970] 
using ISIS I data, Shaw and Gurnett [private communication, 1975], 
and Hosier et al. [1973] using Imp-6 data. 

Hosier et al. [1975] have reported observations of UHR noise 

from frequencies somewhat below f to f.„_» At the time of these 
^ pe UHR 

250 to 600 kHz observations Imp 6 was inside the plasmasphere. Shaw 
and Gurnett [private communication, 1973] have observed the electric 
field of UHR noise at 100 and 178 kHz inside the plasmasphere with 
another experiment on the Imp-6 spacecraft. Their experiment is cali- 
brated to give electric and magnetic field spectral densities. They 
report typical electric field spectral densities of 

^ = 10"^^ to 10"^^ V^/m^z 
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for 100 and 178 kHz noise in the plasmasphere at L values from 4.5 
to 5 close to the equatorial plane, hut the magnetic field of UHR 
noise is apparently too weak to be observed. 

Alouette I and II have also observed UHR noise. The upper- 
frequency cutoff of the noise is generally at or near with the 

lower frequency-cutoff sometimes near max (f , f ) and sometimes 

p0 g0 

near min (f , f ) for both f > f and f > f . The noise is 
pe^ ge' pe ge ge pe 

as large as 50 — 40 db above the background, but the noise level has 
not been reported in absolute terms. Barrington et al. [1971] show 
two samples of the noise intensity as a function of frequency from 
0.2 to 7 MHz near the auroral zone in their Fig\ire 1. Recording 
(a) was obtained near \ = 75° with characteristic frequencies esti- 
mated by Barrington et al. to be f = 0.700 MHz, f = 0.784 MHz, 

P0 g0 

and fyjjp = 1.04 MHz. Noise above background was observed from about 

0.850 MHz to 1.07 MHz. Recording (b) was obtained near \ = 66° where 

Barrington et al. estimated f = 0.182 MHz, f = 0.517 MHz, and 

pe ge 

fyjjR = 0.545 MHz. Noise limits were about 0.15 MHz and 0.658 MHz. 

Cerenkov emission has been suggested as the mechanism respons- 
ible for UHR noise by Walsh [1964], Bauer and Stone [ 1968 ], Gregory 
[ 1969 ], Muldrew [1970], Laaspere et al. [I97l]> Uim and Laaspere 
[ 1972 ], and Mosier et al. [1975]. 

As with VIF hiss, large indices of refraction are required 
for the production of large powers in the UHR noise frequency range 
by low-energy electrons through the Cerenkov mechanism. Large 
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indices of refraction are required for the production of large powers 
in the UHR noise frequency range by low-energy electrons through the 
Cerenkov mechanism. Large indices of refraction exist for wave nor- 
mal angles near the resonance cone angle for frequencies 

^^pe^ ^ge^ < f < fujjjj. Typically there are four or more orders 
of magnitude difference in power production between frequencies in 
this range to those out of this range. The region of expected gener- 
ation of Cerenkov noise from the ambient energetic electrons in the 
magnetosphere, that is, the region for which max ^ge ^ 

for a frequency of 178 kHz, is shown in Figure l8. This region 
corresponds to region 3 of the CMA. diagram [Stix, I962] bounded by 
S = 0 (f = fyjjp), R = (f = fgg), and P = 0 (f = f^^). 

Calculations in the model magnetosphere described above have 
been made at equatorial latitudes in the 100 and 178 kHz generation 
regions to see if the theoretical electric field spectral densities 
are near the ones observed by Imp 6. An energetic electron spectrum 
reported by Schield and Frank [1970] to be representative of the 
plasmasphere and assumed to be isotropic was used in the calculations. 
The spectrum is shown in Figure 3. For a single electron Cerenkov 
waves are generated with k vectors forming a cone around the earth’s 
magnetic field with an opening angle of twice the wave normal angle. 
For a spectrum of electrons a representative wave normal angle is 
chosen, the one corresponding to that part of the spectrum with the 
largest power production. For an isotropic pitch angle distribution 
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of electrons there will he two representative cones, for electrons 
with pitch angles less than and greater than 90°. 

The two-dimensional ray tracing program first determined the 

— f 

ray paths of those waves with k vectors in the meridian plane. The 
two rays which initially began decreasing in altitude went down until 
their wave frequencies f became very close to at >diich point 

they turned around and began increasing in altitude until f became 
close to The attenuation of the waves after the first turn 

was from 2 — 4 db and 10 — 20 db ~ 200 km from the point where 
f = Tyjjp. The two rays which initially began increasing in altitude 
continued until f became close to where the atten\iation again 

began increasing rapidly to over 10 db. Illustrative ray paths in 
two dimensions are shown in Figure 19 for a frequency of IJd kHz and 
a generation point on the equator at an altitude of 7100 km. The 
three-dimensional ray tracing program was then used to calculate ray 
paths for those initial k vectors out of the meridian plane. The ray 
paths for the upgoing rays were very similar to those in the meri- 
dian plane except that they moved in longitude as well as altitude 
and latitude. The initially downgoing nonmeridional rays, however, 
generally followed the same path but reversed directions, or re- 
flected, at higher altitudes. They then continued upward mtil f 
approached where attenuation again began increasing rapidly. 

These ray path calculations mean that at a point in the magnetosphere 
near the equatorial plane frequencies close to should be strong- 


est, decreasing toward f 


pe 
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2 

The volume of space illuminating a one m antenna in the UHR 
2 

region is a one m tube along the ray paths through that antenna 
since the region is fairly •uniformly producing energy. To estimate 
power, electric and magnetic field spectral densities at a point in 
the UHR region, the contributions to the spectral density from points 
along the ray paths through the point are added, including the 
effects of attenuation along the path. 

The fields were calculated at the observer’s position from 
powers generated along the ray paths and added, giving the total 
field at the obser^vation point due to all of the possible generation 
points. Approximate spectral densities at 178 kHz are found to be 

II = 10“^^ w/m^Hz , 

^ = 10-5 _ 

m-52 .„2/ 

=10 Wb /m Hz , 

and at 100 kHz 

dP ,--15 / 2„ 

^ = 10 ^ w/m Hz , 
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^ = 10"^^ Wb%V 

df ' 


The electric field spectral densities are considerably larger than 
those measured by Imp 6. But the magnetic field spectral densities 
are near the noise levels of the magnetic field spectrum analyzers 
on Imp 6 of 

= 2.6 X 10"^^ Wb^/m^Hz 

178 kHz 



and 


= 5.4 X 10"^^ Wb^/m^Hz 

100 kHz 

[R. Shaw, private communication, 1973]. 

There are at least three factors which may explain this appar- 
ent high electric field spectral density. 

(1) The electron spectrum used may have been too Intense which would 
reduce all three i^ectral densities. 

(2) Due to thermal effects or wave normal scattering by irregularities 
some or all of the wave noimial eingles of the waves may have been 
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dispersed away from the resonance cone. Some initial wave nor- 
mal angle dispersion will occur natiirally from the range of vy 
in the energetic electron spectrum. A change of one degree from 
the resonance cone angle will reduce the electric field spectral 
density by about two and one-half orders of magnitude and raise 
the magnetic field spectral density by the same factor. 

(5) Thermal effects would also influence the power production. 
Trulsen [1971] notes that thermal effects dominate near the 
hybrid resonances. A finite temperature will reduce the maxi- 
mum index of refraction, increasing the minimum particle velo- 
city that can radiate. The resonance will also be spread in 
frequency, which may also explain the existence of power outside 
of the limits of the frequency range expected from cold plasma 
theory, max (f , f ) < f < frm-p. But account has not been 
taken of these thermal effects. 

The electric field spectral density was calculated at points 
in the magnetosphere close to where the two recordings mentioned 
earlier [Barrington et al., 1971] were taken by Alouette II. The 
critical frequencies for the first calculation were f = 0.688 MHz, 
fgg = 0.817 MHz, and = 1.07 MHz. The calculated electric field 

spectral density from 0.82 to 1.06 MHz is shown in Figure 20. The 
maximum spectral density of 5-8 X lO" v /m Hz is again considerably 
above the Iirrp-6 measurements. The only number that can be derived 
from Alouette reports is that an average maximxom electric field is 
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4U db above the noise level at 200 kHz. Barrington et al. [1971] 

-8 

give a receiver bandwidth of 4o kHz and a noise level of 7*5 X lO” 
V/m. 44 db then corresponds to a spectral density of 2.5 X 10 
V^/m^Hz. The spectrum shows no obvious similarity to recording (a). 

For the second case the critical frequencies were 
f = 0.177 MHz, f = 0.492 MHz, and f,,^ = 0.523 MHz. Figure 21 
shows the calculated electric field spectral density which is far in 
excess of the measured values. The three possible reasons given 
above for the difference between the calculated and observed spectral 
densities may be valid here too. 

Cerenkov radiation from typical magnetospheric electron spec- 
trums seems more than adequate to provide the power observed in UHR 
noise. The spectrums of calculated UHR noise are not obviously like 
those of observed UHR noise. Thermal effects on the generation and 
propagation of the waves may be responsible for the lower observed 
electric fields. 
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VIII. ELF HISS 


ELF hiss, spectrograms of which are shown in Figure 22, is 
neirrow-band (a few hundred to a few thousand Hz), white noise com- 
monly observed with low-altitude satellites from 50° to 80° invariant 
latitude and 6 to l8 hours magnetic local time [Taylor and Gurnett, 
1968]. Gumett et al. [1971] found that ELF hiss was generally 
observed to be downgoing, but occasionally i;5)going at invariant 
latitudes less than 60°. They found that the transition region 
between upgoing and downgoing ELF hiss was often the boundary between 
the plasmasphere and the light ion trough. They tentatively con- 
cluded that upgoing ELF hiss may be generated at high altitudes out- 
side the plasmapause and reflected upwards at some point below the 
satellite after crossing the plasmapause boundary. 

An intense ELF hiss event from Injun 5 has been analyzed to 
determine its approximate spectral densities. The orbital parameters 
at the time of largest intensity were: altitude = 2174 km, invariant 

latitude = 47-6°, magnetic local time = 15.6 hours. The measured 
wide band field intensities were E = 50 mV and B = 6.5 my. The 
bandwidth was estimated to be 500 Hz and a cosine correlation of 0.5 
was assumed. Spectral densities of 

^ = 1.0 X 10"^ V^/m^z , 



70 


^ = 1.4 X 10"^^ Wb^/m\z , 

^ = 1.5 X 10"^® w/m^z 

were calculated. These spectral densities appear to be the strongest 
spectral densities to have been reported for VLF emission. ELF hiss 
was found to be the strongest VLF emission observed with Inj\in 5 also 
[Taylor and Gumett, I968]. 

The white noise character of ELF hiss suggests that Cerenkov 
emission may be its generation mechanism. However, the facts that 
Cerenkov power is generated near the resonance cone and that down- 
going waves with such wave normal angles reflect when their wave 
frequency approaches fTtm seem to rule out the Cerenkov mechanism. 

The lower- frequency cutoff of ELF hiss was determined by Gtirnett and 
Bums [1968] to be usually below the proton gsrrofrequency, 
which is a considerably lower frequency than fjjjg* The observed 
power spectral density of ELF hiss is also two and one-half orders 
of magnitude above the largest power expected from perfectly guided 
incoherent Cerenkov radiation from an auroral zone flux tube from 
the results shown in Figure 4. In addition, Fig\are 8 shows that the 
largest contribution of 1 kHz power is from eqmtorial regions. Even 
if the power were generated at angles other than those near the 
resonance cone, ray path calculations show that 8OO Hz waves will not 
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propagate to satellite altitudes from initial altitudes greater than 
19^000 km starting on the 65° invariant latitude field lines. At 
15,000 km initial wave normal angles from - 10° ^ ^ 0° will pro- 

pagate to satellite altitudes and this range widens as the initial 
altitude decreases, to - 60° ^ t ^ U0° at an initial altitude of 
5000 km. 

It was thought that ray path calculaticns of waves with ELF 
hiss frequencies would be able to confirm the theory of Gumett et al. 
[1971] predicting the reflection of downgoing ELF hiss at low alti- 
tudes near the plasmapause. They predict that ELF hiss waves would 
travel to significantly lower latitudes while reflecting at low 
altitudes, resulting in the occasionally observed change from down- 
going outside the plasmapause to upgoing within the plasmapause. 

The two-dimensional ray tracing program was used to calculate ray 
paths for initial altitudes above satellite altitudes and for invar- 
iant latitudes slightly larger than 60° (the invariant latitude of 
the L = 4 plasmapause ) . The rays did reflect at the L = 0 cutoff as 
predicted by Gumett and Burns [1968] but the decrease in latitude 
was never more than a few degrees. Gumett et al. [1971] show an 
exang)le of ELF hiss which changes direction (their Figure 5) at an 
invariant latitude of about 59° • Inspection of the original data 
shows that the ELF hiss was upgoing at invariant latitudes from 38° 
to 59° and downgoing from 59° to 69°. The satellite altitude and 
magnetic local time stayed nearly constant during the ELF hiss event 
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at about 2400 km and 6 , respectively. For the ELF hiss waves to be 
observed downgoing and then upgoing at least some waves in this event 
would have had to propagate 20° in invariant latitude (~ 1000 km) 
at altitudes below 2400 km. The magnetosphere model used here will 
not reproduce such behavior, but perhaps occasionally there are hori- 
zontal ducts to produce large latitudinal displacements of ELF hiss 
propagation paths. 

An alternative explanation of the behavior may be that the 
generation region is below satellite altitudes in the plasmasphere 
and that the waves propage upwards and then reflect at high altitudes, 
coming down outside the plasmasphere at higher latitudes. Two objec- 
tions to this theory are immediate, however. First, the generation 
region would necessarily be small implying a very large volume 
emissivity in a region without large electron fluxes. Second, there 
appear to be no plasma cutoffs or resonances^ capable of reflecting 
upgoing waves above satellite at ELF frequencies. Ray path, calcula- 
tions confirmed that there are no reflections of upgoing ELF hiss 
frequency waves starting in the plasmasphere above satellite alti- 
tudes. 

On the basis of the observations ^d the calculations pre- 
sented here, it seems most likely that the generation mechanism for 
ELF hiss is a particle cooperative one, for exan5>le coherent Cerenkov 
radiation or a wave-particle instability. The generation region 
appears to be at altitudes above a few thousand kilometers but below 



75 


a few tens of thousands of kilometers near the plasmapause or at 
slightly higher latitudes. The waves propagate down the field lines 
and reflect at the L = 0 cutoff, sometimes propagating to lower 
Latitudes at low altitudes. 
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IX. DISCUSSION 


A. Applicability of the Theory 
There are two assumptions in the Cerenkov theory which should 
be closely examined, the assuptions that the plasma is cold and that 
the test particles and the radiation they emit have no effect on the 
plasma. Trulsen [1971] states that thermal effects are particularly 
important in a plasma near the hybrid resonances where = 0 in cold 
plasma theory, in 5 )lylng that thermal effects are important in the 
consideration of VLF hiss and UHR noise. A finite temperature may 
be introduced into the theory by using the Vlasov equaticn to deter- 
mine the dielectric tensor describing the plasma [Trulsen, 1971] or 
by introducing a temperat\ire dependent pressure term in the equation 
of motion. In either development does not go to zero at the 
hybrid resonances. Because of this fact Trulsen [l97l] claims that 
the Cerenkov power spectra is changed ’completely’ from the cold 
plasma result near fTTm^ removing the singularity in the spectra at 
f = and extending the spectrum far above ^ finite tempera- 
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ture will, in smearing out the resonance, reduce the maximum index of 
refraction n. This reduction in the maximum n will reduce the power 
emitted from an electron spectrum by increasing the minimum parallel 
particle velocity that will radiate. 

Some theoretical calculations may be made to determine the 
necessity of including thermal effects in the power and ray tracing 
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theories. For cold plasma theory to be valid 3 (gas energy density/ 

magnetic energy density) should be much less than one. In an extreme 

case magnet ospheric plasma parameters might be = O.h G, 

n = 10^ cm”^, and T = 10^ which gives a P of 2 X lO"^. !^e com- 

-¥ 

ponent of the phase velocity of the emitted wave along v^^^ 

should be considerably higher than the mean electron thermal velocity. 

With the choice of collisions as a dissipative force and the numbers 

used in the calculation, the thermal velocity was as large as v . 

' Ph> II 

Perhaps some of the discrepancies that do exist between theory and 
observation presented here are explicable in terms of thermal effects. 
Because of the much more complicated dispersion relation it was 
necessary to assume zero temperature in order to obtain a reasonable 
numerical description of the plasma, however. 

The assumption that the particles and the emitted waves have 
no effect on the plasma may also be considered. It will be shown 
later in this section that the conversion efficiency from particle 

_7 

energy to VLF hiss wave energy is no larger than about 10 . Therefore 
the wave energy is much less than the particle energy, so that neither 
the wave emission nor the propagating waves can affect the particle 
spectrum. In the case of VLF hiss it may be that the emitted waves 
do interact with the plasm and energetic particles, causing coherent 
emission which greatly increases the wave energy output. (This 
possibility will be briefly considered later.) 
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B. Validity of the Programs 


The program which calculates the power emitted by a particle 

has been checked against published data of, similar calculations by 

Jorgensen [1967] and Lim and Laaspere [1972 ]• In these checks the 

collision frequency was taken to be zero to agree with the previously 

reported calculations. J(^rgensen's Figure 10 was the basis of the 

first coniparison. With the plasma parameters given for an altitude 

4 

of 10 km, the power was calculated for five frequencies and was 
within about 5^ of the value read from the figure. 

Figure 1 of Lim and Laaspere [1972] provided the second com- 
parison for the power program. Again using the plasma parameters 
listed, at a frequency of I/2 f^^ the power was calculated for four 
energies. The calculations and the data from Lim and Laaspere 's 
figure agree to within about 5^- 

Three tests were made of the ray tracing program to insure its 
validity. The first test was for a ray path in a horizontally 
stratified medi\im. The plasma density and magnetic field were 
assumed to vary only in the vertical direction, implying that the 
index of refraction varies only vertically. A coordinate trans- 
formation from rectangular to polar coordinates was used to make full 




78 


degrees from the resonance cone angle. At this point rovindoff errors 
are expected to be large. Varying both the magnetic field and the 
plasma density resulted in constant to within 0.008^ for -fre- 
quencies more than 0.04^ away from f^^, where roundoff errors are 
also expected to be large. 

It was shovm by Shawhan C1966] that if a ray is reversed it 
must retrace the same ray path. This characteristic of ray paths 
provides a means of checking the integration accuracy of the program. 
This reversibility check was made, and for a two-way distance of 
16,000 km, the ray returned to within 0.001^ of its starting point, 
compared to the one-way distance of 8 OOO km and to within 0.001^ of 
its initial wave normal angle. 

Finally, the ray tracing programs were compared to that of 
Kimura [I966] by reproducing the ray path he calculated in his Figure 
6. The plasma density model involved was an electron- three-ion 
diffusive equilibrium model and the magnetic field was represented 
by a centered dipole. For a frequency of 1 kHz, an initial altitude 
of 500 km, an initial latitude of 30®, and an initial wave normal 
angle (with respect to the radius vector) of 0®, ray paths were cal- 
culated. The turning points (latitude maximums) found by Kimura were 
given in his Table 4. The turning points found here agree to within 
a few percent in altitude and latitude of the ones fo\uid by Kimura. 
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C. Coherent and Amplified Cerenkov Radiation 
The Cerenkov power calculation presented here for several 
magnetospheric wave phenomena, as in Jorgensen [I 968 ] and Lim and 
Laaspere [1972] for VLF hiss, assumed that each electron radiated 
independently, that is, with random phase. The total power rad- 
iated incoherently by n particles is 

Pt = nP 

where P is the power radiated by each particle. On the other hand, 
if the n particles radiated coherently, that is, all with the same 
phase, the total power would be 



Therefore, if the electrons are radiating somewhat coherently in the 
Cerenkov mode an approximately square law relation between electron 
and VU* fluxes would be expected since the radiating power is pro- 
portional to the square of the number of radiating electrons. It 
is clear that a relatively small number of electrons emitting Cerenkov 
radiation coherently could substantially increase the power to the 
observed levels. 

The threshold effect in the relationship between low-energy 
(~ 100 ev) electrons and VLF hiss fomd by Gumett and Frank [1972] 
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for VLF hiss might he due to the steepness of the square law being 

mistaken for a threshold or it may be that there is a threshold cor- 

. U 5 / 2 

responding to the onset of an instability at 10 —10 elec/cm sec ev 

ster. The observation by Gurnett and Frank [1972] showed VLF hiss 

power changing about 3 l/2 orders of magnitude while the particle flux 

changed by about 1 l/2 orders of magnitude \diich is approximately 

square law, which suggests partially coherent Cerenkov radiation. 

Singh [1972] has foxind that Cerenkov waves may be an?)lified by 
a factor up to three orders of magnitude by interaction with an elec- 
tron beam. Significant amplification factors for an interaction 
length of 1000 km were found in the frequency range 0.1 f^ ^ f f^, 
where f^ is the plasma frequency of the beam. For a beam density 
of 1 cm ^ (typical for the spectra in Figure 5), Cerenkov waves would 
be amplified in the frequency range 0.9 ^ f ^ 9 kHz. Singh's calcu- 
lations, however, assumed a longitudinal dispersion relation. If 
similar results are valid for wave normal angles near the resonance 
cone angle, VLF hiss may be explicable in terms of amplified Cerenkov 
radiation. 


D. Instabilities 

There are a nximber of plasma instability mechanisms that have 
been proposed to explain various magnetospheric wave phenomena. 

Kennel and Petschek [1966] have determined that the whistler 
mode is unstable' to wave growth when the energetic electron pitch 
eingle distribution is sufficiently anisotropic, favoring perpendicular 
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velocities. The whistler mode noise is generated at high altitudes 
near the equatorial plane at the electron gyrofrequency and may he 
Doppler shifted into the VLF frequency range. 

The generation mechanism for VLF and ELF hiss and other pheno- 
mena may involve generation of waves in the drift mode at frequencies 
between the ion and electron gyrofrequencies and subsequent coupling 
to the whistler mode [Sizonenko and Stepanov, 1967] • The fastest 
growing whistler mode has a frequency given by 


u 


f = 




ge gi 


where u^ is the ion velocity perpendicular to due to the influence 

of an ion cyclotron wave and v^ is the ion thermal velocity. In the 

X/2 

high plasma frequency limit (f f.) is the lower hybrid reso- 

nance frequency (f,.— ) (which is usually assumed to be the lower- 

IjHK 

frequency cutoff of VLF hiss). Therefore this fastest-growing 

whistler mode is in the correct frequency range for VLF phenomena. 

Lee, Kennel, and Kindel [1971] found a high-frequency (near 

f,._) Hall current instability due to electron and ion drifts perpen- 
ItlK 

dicular to B^. They determined that this instability could produce 
waves in the ELF and part of the VLF frequency band. 


Further consideration of these generation mechanisms is not 
undertaken here, but it seems clear that many of the wave phenomena 
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must be due to cooperative plasma effects, 
the instabilities mentioned above. 


perhaps one or more of 
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X. JUPITER 

Jupiter is known to emit intense bursts of decametric radia- 
tion up to 10 ^ w/m Hz for frequencies less than 40 MHz [Carr and 
Gulkis, 1969 and Warwick, 1967 ]. Some of this radiation is associated 
with the orbital position of the satellite lo and is made up of micro- 
second bursts which drift in frequency. This frequency and time 
structure might be attributed to particles passing through a region 
with varying plasma parameters and emitting radiation characteristic 
of the plasma [Warwick, I 965 ]. 

Assuming that the characteristic frequency is the electron 
gyrofrequency, incoherent Cerenkov (s = O) and cyclotron (s = 1, 2) 
power calculations were made at 26 MHz. A model Jovian magnetosphere 
with a centered dipole magnetic field was assumed to give a maximum 
gSrrofrequency at the latitude of lo's field line at the Jovian sur- 
face of about 40 MHz and the magnetosphere was assumed to be a 
hydrogen plasma in diffusive equilibrium. An estimate of the ener- 
getic electron spectrum emitted and accelerated by lo [R. Hubbard, 

1972 ] up to 300 keV was used in the power calculation. Estimating 

_2 

the r power loss from the initial spread in ray angles the received 

-35 

power at the surface of the earth was calculated to be about 10 

/ 2 l4 

w/m Hz. This is a factor of about 10 less than observed during 


decametric bvirsts. The power from the s = 0 and 1 modes were nearly 
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equal and the power from the s = 2 mode was a factor of two lower 
than the power from the s = 1 mode. The observed rapid tinie varia- 
tions and this power deficit calculated for incoherent radiation 
sijggest that the emission mechanism is a coherent one. The coherence 
in time or space might be expected to be short, which would give a 
bursty character to the emitted radiation. 
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XI. CONCLUSIONS 

The calculations of the power of incoherent Cerenkov radiation 
and the ray paths of resulting waves have been helpful in vinder- 
standing the characteristics of some broadband LF, VLF, and ELF ob- 
servations near the earth. Qualitative and quantitative comparisons 
between the theory and observations of VLF hiss, saucers, fast 
hisslers, UHR noise, and ELF hiss have been made. Calculation of 
decametric radiation from Jupiter have also been made. The results 
of these calculations are summarized below: 

(l) VLF hiss. A calculation of expected VLF hiss power using an 
electron spectrum measured by Injun 5 was compared to a simul- 
taneous observation of VLF hiss on the same spacecraft. The cal- 
culation included collisions and an estimate of thermal effects. 
Ray path calculations from the presumed generation region above 
the satellite were used to detennine theoretical spectral densi- 
ties and shapes. The power and magnetic field calculated were 
below observed levels by factors of 10^'^ and 10^, but the elec- 
tric field calculated a factor of three above observed levels. 

From Injun 5 data, VLF hiss waves were determined to be 
propagating at wave normal angles about six degrees from the 
resonance cone. When the spectral densities of the fields were 
calculated for such a wave normal angle they were each about two 
orders of magnitude below the observed spectral densities. The 
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calculated spectral shape of VLF hiss agreed well with observed 
spectra, however. This agreement was due to larger latitudinal 
and longitudinal spreading of the higher frequency waves from a 
source region confined to the flux tube of the electron precip- 
itation region. Nevertheless, it seems unlikely that VLF hiss 
is generated by incoherent Cerenkov radiation from electrons. 

The calculated spectral shape of V-shaped VLF hiss is 
similar to observations when a source region limited in altitude 
extent at approximately 4000 km altitude was assumed. However, 
because of the small emitting volume the output power was low, 
which reduces the likelihood that the observed power is generated 
by completely incoherent Cerenkov radiation. 

( 2 ) Saucers. Ray path calculations defining the spectral shape of 
saucers tend to confirm the theory that the generation region is 
below the satellite and that the emission is caused by upcoming 
electrons. The characteristic size of the source region must be 
only a few hundred kilometers since the envelope of the saucer 
is well defined. The generation mechanism appears not to be in- 
coherent Cerenkov due to the low power available from the small 
source region. The initial wave normal emgle for the emission 
mechanism appear to be very close to the resonance cone angle. 

(3) Fast hisslers. Ray path calculations show that the generation 
regions of fast hisslers are at high altitudes, e.g. 15,000 km. 
Incoherent Cerenkov radiation produces enough power with a gen- 
eration region size of a few thousand kilometers to match the 
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observations of fast hisslers, but the short emission time and 
the apparent inability of Cerenkov waves (with initial wave nor- 
mal angles near the resonance cone) to reach the ground from 
these altitudes argue against the Cerenkov mechanism. Only a 
small range of initial wave normal angles nearly aligned with 
the magnetic field (0.2° — 2.3°) initially will reach the ground 
from an altitude of 15,000 km. It might be that this Cerenkov 
noise is momentarily scattered into a downgoing wave normal 
range, but this irregular scattering structure has not been 
considered. 

(4) UHR noise. Observations of UHR noise show wider spectral shapes 
and lower spectral densities than calculated from incoherent 
Cerenkov radiation. Thermal effects may reduce the spectral 
density and change the spectral shape by increasing the minimum 
energy vdiich an electron must have to Cerenkov radiate and by 
altering the wave normal angle distribution. Thermal effects 
would broaden the resonance, reducing the maximum possible index 
of refraction and increasing the range of emitted wave normal 
angles. This change in the dispersion relation would effect pro- 
pagation as well as generation. Since the calculated power was 
above the observed power, it is possible that a calculation 
including thermal effects would predict well the characteristics 
of UHR noise. 

(5) ELF hiss. Observed power levels and wide ranges of occurrence 
in latitude and longitude for ELF hiss are not consistent with 
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calculations of incoherent Cerenkov radiation. In addition, for 
these frequencies wave normal angles near the direction of the 

3 

magnetic field are required for wave propagation from hi^ alti- 
tudes to satellite altitudes. 

(6) Jupiter. Incoherent Cerenkov (s = O) and cyclotron (s = 1, 2) 

I 

radiation from electrons emitted from lo and accelerated hy lo's 
sheath electric field was found to be fourteen orders of magni- 
tude less than observed during decametric- bursts. Most of the 
power was calculated to be propagating perpendicular to Jupiter ' s 


magnetic field. 



Observed and Calculated VLF Hiss Intensities 
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Figure 1 


Figure 2 


Figure 5 


Figure U 


Figure 5 


Figure 6 


Figure 7 


Figure 8 


Three representations of an index of refraction surface 
^LHH ^ f ^pe- 

A spectrogram showing VLF hiss when the satellite was near 
an auroral arc [after Mosier and Gurnett, 1972]. 

Energetic electron energy spectra used in previous and 
present calculations of magnet ospheric noise generated 
by the Cerenkov mechanism. 

The spectral density of observed VLF hiss and of calculated 
incoherent Cerenkov radiation by various investigators. 

The variation of wave normal angles over a ray path as a 
function of altitude for downward propagating VLF hiss 
waves. The wave normal angles stay very close to the 
resonance cone angle over the entire path. 

The latitude dispersion of VLF hiss ray paths for frequencies 
of 2 and 10 kHz and two generation altitudes on the 70° in- 
variant latitude field line. 

The electron density derived from an electron-three ion 
diffusive equilibrium model of the magnetosphere with a 
plasmapause at L = 4. 

The incoherent Cerenkov power emitted in altitude segments 
up to 48,000 km altitude at 70° invariant latitude. 
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Fig\are 9 


Figure 10 


Figure 11 


Figure 12 


Figure 15 


Figure l4 


Power and electric and magnetic field spectral densities 
expected from Cerenkov radiation ft*om auroral electrons at 
an altitude of 2500 km. 

The dependence of electric and magnetic field spectral 
densities upon wave normal angle for a typical VLF hiss 
wave at 2500 km. 

The emission intensity for an auroral electron spectirum for 
various modes. The s = 0 (the Cerenkov mode) dominates 
from f_^ to frequencies slightly below f . The s > 0 

LiriK p© 

modes are called normal cyclotron modes and the s < 0 modes 
are called anamolous cyclotron modes. 

The spectrogram depicting downgoing V-shaped VLF hiss from 
22h 54m 50s to 22h 57m 50s and an upgoing saucer centered 
at 22h 58 m 00s after Gurnett and Frank [1972] (originally 
in color in their Plate 5)* 

The latitudes at which theoretical V-shaped VLF hiss waves 
pass an altitude of 2500 km for a generation region at an 
altitude of 4000 km and invariant latitude of 73°. 

i 

A spectrogram showing an intense upgoing saucer after 
Gumett et al. [1971] (originally in color in their Figure 


6 ). 


Figure 15 The latitudes at ■vrtiich theoretical saucer waves pass an 

altitude of 2500 km for a generation region at an altitude 
of 1800 km and an invariant latitude of 69 * 5° • 
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Figiire l6 


Figure 17 

Figure l8 
Figure 19 


Figoire 20 
Figure 21 
Figure 22 


The range of initial wave normal angles that will propagate 
into the earth-ionosphere waveguide for a generation alti- 
tude appropriate for fast hisslers (l5,000 km) at an im'ar- 
iant latitude of 71° (Byrd Station). 

Theoretical and experimental fast hissler dispersions from 
2—12 kHz. 

The 178 kHz cold plasma UHR noise generation region. 

178 kHz UHR noise ray paths from a isotropic electron spec- 
trum at an altitude of 7100 km on the geomagnetic equator. 

Note the reflections at f and the attenuation of the 

pe 

waves. 

A theoretical UHR noise electric field spectrum near 1 MHz 
in the aiiroral zone. 

A theoretical UHR noise electric field spectrum near 0.5 
MHz in the a\iroral zone. 

A spectrogram of ELF hiss observed with the Injun 3 satel- 


lite. 



INDEX OF REFRACTION SURFACE FOR F = 8 KHz, ALTITUDE = 2500 KM. F^^r = 5.99 KHz 
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Figure 7 




















C-G73-273-* 


121 



Figiire 10 








o 

m 




> 

Id 

tr 


^ ^^ 12 $ 




I . t 


AVx I 

V ') ^^> />r- ^fy V|^ 


uS ^ 


CM CO 


CM 

O 


CM CD CD O 


CM 


03 rvi ^ 

^-) — CNJ to 

C\i K CD <J> 
C\i <i> — — 


<x> „ ^n 

to C£» — ; Kf 
CM — CD O 
CM D- — CM 




I — ^ 

3 X liC 

^ m m 

to CO cn a:) 

CM fO CO O 
CM N- . — CM 


> h h 

X -M -M 

- 2 < 


O 'O o 

{ZH>i) A3N3n03fc!j 


Figure 
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Figure 15 



FREQUENCY-kHz 
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Figure l6 
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Figure 17 
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Figure l8 
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ALTirUDE-km 


Figoxre 19 
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Figure 20 
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Figure 22 



